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stations in which the cable tension is zero (non- spinning case) are not 
treated herein. Differential equations of vibration of the elastic 
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motions are coupled through the nonhomogeneous boundary conditions of 
the cable. This mathematical difficulty is resolved by using the con- 
cept of concentrated fictitious masses. The cable equation is solved 
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DEFINITION OF SYMBOLS 


Symbol 

c 

El 

I 1 y I2>l3 

£ 

L 

m 


M*(t) 


P(t), s(t) 
ri(x) 

R 

t 

u(x, t) 
v,w 


x,y,z 

X, Y,Z 

XI, Xp,X 3 

Pn L 

\ 

\x 

p 


Definition 

center of rotation to space station 
bending rigidity of cable 

principal moment of inertias of space station 

distance of cable connection to c.g. of space station 

length of cable 

mass of space station 

applied moment 

generalized coord inates 

eigen-function of a free-free beam 

center of rotation to c.g. of space station 

time variable 

displacement relative to moving coordinates 

components of u relative to moving coordinates 

rotating coordinates 

inertial coordinates 

body- fixed principal axes 

constant in eigen-function of beam 

eigen-value of free- free beam 

eigen-value of matrix (A -1 !!) 

eigen-value of matrix Q 

mass density per unit length of cable 


DEFINITION OF SYMBOLS (Continued) 


Symbol Definition 

^,9,0 Eulerian angles of space station 

ty, 9,0 column matrices of the Eulerian angles 

Barred symbols refer to counter-weight of space station. 

Subscripts i, n denote the ith and nth mode of vibrations of cable. 
The dot above a symbol denotes the time derivative. 
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TECHNICAL MEMORANDUM X- 53650 


ON DYNAMICS OF TWO CABLE- CONNECTED SPACE STATIONS 

SUMMARY 


Presented is an analysis of the dynamic problems of two cable- 
connected space stations rotating about an axis normal to their orbital 
plane to provide artificial gravity. The dynamics of cable-connected 
stations in which the cable tension is zero (non- spinning case) are not 
treated herein. Differential equations of vibration of the elastic 
cable and the angular movements of the stations are derived. These 
motions are coupled through the nonhomogeneous boundary conditions of 
the cable. This mathematical difficulty is resolved by using the con- 
cept of concentrated fictitious masses. The cable equation is solved 
by using Galerkin 1 s approach for both free and forced oscillations. A 
general nth order de terminantal frequency equation of free vibration of 
the system is obtained. The responses of the space stations to applied 
time-varying moments are presented in analytical forms. 


I. INTRODUCTION 

The dynamic problems of cable-connected space stations rotating in 
the orbital plane have been discussed in much of the published literature. 
Reference 1 presents a detailed analysis of many types of arrangements 
of the stations, as well as a bibliography on this subject. The con- 
figuration of the space stations treated in this report is two arbitrarily 
shaped space vehicles which are connected by a long elastic cable. The 
system is considered rotating about an axis normal to the orbital plane 
with constant angular velocity. The dynamic problems concerned here are 
(1) the dynamic stability criteria, (2) the natural frequencies of free 
vibrations and (3) the dynamic responses of the stations to applied 
moments . 

In a near-zero earth-gravitational field, the dynamic behavior of the 
space stations in free flight is predominated by the response of the con- 
necting cable. Therefore, it is essential to develop the dynamic equation 
of the rotating cable coupled with the free angular movements of the space 
stations. This concept leads to an approach entirely different from that 
used in reference 1. 



The equations of motions of the space stations and of the elastic 
cable are formulated independently, but the two sets of variables are 
linked kinematically at the points of connection. This enables us to 
eliminate the angular variables of the stations. By using concentrated 
fictitious masses to represent the dynamic coupling, the problem is 
reduced to the solution of the cable equation with nonhomogene ous mass 
distribution but homogeneous boundary conditions. 


II. ANALYSIS 


A. Description of the Coordinate Systems 

The formulation of the dynamic problems of space vehicles can be 
greatly simplified by carefully choosing the coordinates and the vari- 
ables. To achieve this, several coordinate systems are required (see 
figure 1). 

(1) The inertial coordinates XYZ: Let the X and Y axes be in 
the orbital plane of the space stations and the center of mass of the 
system be the origin. This coordinate system has fixed orientation with 
respect to earth but with moving origin. 

(2) The rotating axes x*yz and xyz: The z-axis is coincidental 

with Z about which the coordinates rotate with a constant angular veloc- 
ity ft. Another set of rotating axes, xyz, with their origin at one end 
of the cable will be used in dealing with the cable equations. 

(3) The principal axes XjX^x 3 and the Euler angles: Let the 

center of mass of the space station be the origin of the principal axes 
and let the x 3 -axis of the undisturbed station be parallel to Z. The 
cable is attached to the space station at a distance jo on x^axis. The 
orientation of the disturbed station can be completely defined by Euler* s 
angles. First, as shown in figure lb, the station rotates about the 
axis normal to the orbital plane with angular displacement \|/. This 
moves the other two principal axes to x^ and x^. Since the two stations 
are rotating as a unit about the Z-axis with angular velocity ft, the 
first Euler angle is ijr + ftt. The second rotation, which is about the 
x'p-axis with angle 0, shifts x^ to its final direction x 1 and the other 
axis to x 3 . The last rotation <j> about x-^axis gives the final orienta- 
tions x 2 and x 3 . 

(4) Symbols for the counterweight; The barred symbols are used 
for the counterpart of the space station which will be referred to as the 
counterwe ight . 
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B. Formulation of the Moment Equations of the Space Stations 


Let us consider that the masses of the space stations are much 
greater than the mass of the connecting cable and that the system is 
free of external forces. Also, we may restrict the disturbed motions 
of the system as follows: 

(1) The motion of the space stations is in the plane of 
rotation and along the line connecting the two c.g.'s, 
and 

(2) the orbital motion of the center of mass of the system is 
not affected by the disturbance. 

We shall now derive the equations of motion of the space stations and 
the connecting cable separately. 

Under the assumption that the cable cannot be extended, the radial 
velocity of the space station is small in comparison with its absolute 
angular velocities co l5 oo 2 and oo 3 about the principal axes. Hence, the 
rotational kinetic energy can be used as the kinetic energy of the 
body: 


T — 2 “I" I3OO3) , (1) 

where I l5 I 2 and I 3 are the principal moments of inertia. Upon sub- 
stitution of the relationships of the absolute angular velocities with 
the Euler angles. 


oox = $> - (Q. + \jr) sin 0 

oo 2 = 0 cos 0 + (H + {jO cos 0 sin $ 

• • 

w 3 = (ft + \|r) cos 0 cos $ - 0 sin <J>, 


( 2 ) 


we have 

T = j Ii[i> - (Q + ijr) sin 0] 2 + j I 2 [0 cos $ + (ft + i) cos 0 sin $] 2 



I 3 [(ft + \|r) cos 0 cos <J> - 0 sin $] 2 . 


( 3 ) 
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Applying the Lagrangian equation of motion, 


_d_ St_ _ 5v 

d t Sq . Sq . . i 9 

i n i n i 


results in, after neglecting the nonlinear terms 7 three moment equations 

m + - ™ + f 


m = 1 2 e + (i 3 - i x ) n 2 e + (i x + ip - i 3 ) n$ + 


dV 

c)0 


(4a) 


M = I x 4» + (I 3 - I 2 ) ff 2 * - (I x + I 2 - I 3 ) no + 


• . c)V 




where M is a moment due to all external forces about the axis through 
the center of mass. Three identical equations can be written for the 
counterweight: 


Mt = I 3 \jr + — 

= i 2 e + (i 3 - n 2 o + (h + - i 3 ) (4b) 

oO 

m- = + (i 3 - i 2 ) n 2 o - (ii + i 2 - i 3 ) no + ~ . 


The potential energies due to the artificial gravity for the two bodies 
are 


V = mRn 2 (-AR) and V = mRQ 2 (-z^) , 


(5) 



these terms become cricical only when the spinning rate is very small. 
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where 


c 


2 


■ZR = — 


) dx + £(l - cos 0) + £(l - cos \[r) 


-ZR = 


) dx + ^ ( 1 - cos 0) + £(l - cos Mr) 

ox 1 


in which u = u(x,t) is the deflection of the cable relative to the rotating 
axes. Notice that the potential energy of the cable is small and can be 
disregarded . 


C. The Vibratory Motions of the Connecting Cable 

For the convenience of using the tabulated data and the integration 
formulas of eigen-f unctions given by references 4 and 5, the origin of the 
rotating coordinates is transferred to one end of the cable, and the new 
coordinates are denoted by xyz and their unit vectors by i,j and k. Let 
M(x,t) be the bending moment due to centrifugal load, N(x,t) the axial 
force in the cable, and \x the position vector of an element from the 
center of rotation, 


u*(x,t) = (c - x) i + v(x,t)j + w (x , t)k . 


The equation of vibration can be expressed by either 


EID 4 u + D^Cx, t) + p(D^ t) = 0 
xx c a 


(6a) 


or 


EID 4 u - D (ND u) + p(Df; t) = 0, (6b) 

x xx t a 

where D x and D t denote partial differential operators with respect to x and 
t, respectively. The last term in equation (6) denotes the absolute acceler- 
ation excluding the acceleration of the origin. By the principle of differ- 
entiation of a vector. 
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( 7 ) 


(Pf <?) = D 2 u + (fc) x ($c) X u + 2(fik) X (D u*) + (tk) X u 

l a. l l 


= - [2nD v + (c - x) n 2 ]i + (D?V - Q 2 v)j + Dfwk. 

t L t 


The formulation of the middle term of equation (6) is given in appendix A. 
Now, the following equations of motions are obtained: 


EID^v + [\ p(x 2 - 2cx) - mR] D 2 !) 2 ^ + p fi 2 (x - c) D v + p(D?v - £2^) = 0 
x z x x t 

(8a) 

EID^w + [1 p(x 2 - 2cx) - mR] Q 2 !) 2 ™ + pfl 2 (x - c) D _w + pDfw = 0. (8b) 

X z. X X t 

There is no known exact solution to the above equations. However, 
an approximate solution can be obtained by using the well-known Galerkin f s 
approach. Fictitious masses are used to provide conditions in which it 
is applicable. 


D. The Boundary Conditions of the Cable 


1 . The Kinematic Boundary Conditions 

It can be seen from figure lb that the angular displacements of 
the space stations and the displacements of the cable are related as 
follows : 


v(o,t) = mt) 

-v(L,t) = I\j/(t) 


w(0, t) = f 0 (t) 
w(L,t) = i0(t). 


(9) 


As a result of these important relationships, it enables us to combine 
the three sets of equations into a single set. 
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2. The Dynamic Couplings 


The external moments in equation (4) come from two sources. 

One source is the applied moments from the control rockets on board the 
space stations as denoted by M'% and the other is the elastic shearing 
forces at the ends of the cable. Conversely, the disturbed motions of 
the stations act upon the ends of the cable as shearing forces. The 
action of the stations on the cable can be considered as concentrated 
masses attached to the ends of the cable with magnitudes satisfying the 
following dynamic conditions: 


\ - * m v <D f v<0 - t »a +M t * » v (Dfv(L, t » a ' M | ( 10a > 


M 0 - /:m w < D jM0,t)) a + m* m; - S w (D^(L,t) a + m| , <iob) 


where m v and m w are the fictitious masses with respect to the motions of 
v and w. Let us assume that the ends of the cable are free to rotate with 
respect to the stations; thus, 


M = M* and M- = . 

0 0 0 $ 


(10c) 


By substituting the right side of equation (4) into equation (10) and making 
use of equation (9), we now reduce the problem to free vibrations of a 
rotating cable with fictitious concentrated masses attached to the ends. 


E. The Modified Eigen-Functions 

Let us consider that the fictitious masses are distributed uniformly 
over a small interval e at both ends of the cable. Consequently, the 
constant p of the last term in equation (8) should be substituted by the 
mass density functions ^(GjX) and p^e^x) with respect to the motions 
v and w. These functions are defined as 


P v (e> x ) 


0 < x < e 

p + m v /e 
P + 


e < x < L - e 
P 
P 


L - e < x < L 

p + m v /e 

p + m I £ 

H w 
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We now represent the two equations of equation (8) by an equation of a 
single variable u(x,t) in the form 


EID 4 u + mfi 2 


P f X p V 

- <2 X cx) - R 


D'-u + - (x - c) D u \ + p (g , x) (D^u) = 0. 

x m N x | u v, ' v ta 


( 11 ) 


1 . Eigen-Functions of a Uniform Free-Free Cable 


The eigen-functions of a uniform, free-free cable which satisfy 
the differential equation, 


d 4 r . 


- - Bfr. = 0, i = 1,2 ... 


dx^ 1 i i 


and the orthogonality condition 


or (x) r.(x) dx = nLS . , (c . is Kronecker delta) 
K n i f ni 5 ni ' 


are 


r i( x ) “ 1 rigid body mode, symmetric 

r 3 (x) = d(c - x) rigid-body mode, anti-symmetric about x ~ c 

d = (c 2 - cL 4- L 2 / 3)” 1 / 2 


r . (x) = cosh 6.x 4- cos 6.x - a.fsinh 6.x + sin 6.x), i ^ 3. 

l 1 i i v 1 l ' l 


Notice that the ith mode given above is the (i - 2)th mode of references 
4 and 5. 

2. The Modified Eigen-Functions Adapted to Equation (11) 

Let r^(c,x) be the eigen-f unction of a free-free cable with a 
mass density function p u (e,x) of which c serves as a small parameter. 

In the limiting case, we have 


lim r . ( g,k) = r . (x), i = 1,2, ... . 

e — » 0 
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Furthermore, the orthogonality condition for the modified eigen-func- 
tions is 


p u (c,x) r n (e,x) r i (c,x) dx = 0, n ± i. 


(12a) 


In dealing with the integration 


P u (c,x) r n (e,x) f(x) dx = / pr n (e,x) f(x) dx 


+ / ( m f /e) r (hx) f(x) dx + / (m /c) r (e,x) f (x) dx, 


L-e 


we may apply the mean-value theorem to the last two integrals of the 
above and obtain 


L 


L 


Lim J p^(c,x) r n (e,x) f(x) dx = J pr^(x) f(x) dx + rn^r^(O) f(0) 
o o 


4- m r (L) f (L) . 
u n N 7 v 7 


Notice that 


(12b) 


r i(0) = 1 r 3 _(L) = 1 

r 2 (0) = dc r 2 (L) = -d (L - c) 

r n (0) = 2 r n (L) = -(-l)" 2 n = 3,4, ... . 
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F. Approximate Solution by Galerki^s Approach 


The general approach to the approximate solution of equation (8) is 
to assume the solution in the form 

00 

v(x,t) = r j .(e,x) p^t) 

i=l 

(13) 

CO 

w(x,t) = r.(e,x) s^t). 

i=l 


In these expressions, r^(e,x) is the modified eigen-function which satis- 
fies the given boundary conditions; p^(t) and s^(t) are the variables to 
be determined* By the same token, we may write 


* (t) = ^ t.(t) 
i=l 


0 (t) = 2 J 0 i (t) 

i=l • 


®(t) = ) %(t) 

i=l 


5r(t) = ) | t ± (t) 

i=l 


e(t) = ) e t (t) 

i=l 


kt) 

i=l 


(14) 


In accordance with equation (9), we have the following kinematic coupl- 
ings between the variables: 


Rr i (0) Pi (t) = £f L (t) 


Rr.(O) Si (t) = J g0 i (t) 


Rr ^ (L) p ± (t) 


-^(t) Rr.(L) s i (t) = ie.(t). 


(15) 
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By making use of equations (4) and (13) to (15), the dynamic boundary 
conditions given by equation (10) become 


iRm r . (0) (p.) 
v r l a 


(I 3 R/i) r (0) p + 

i 


5v 



(16a) 


iRm r . (0) (s.) = (I 2 R/.C) r (0) s + (I 3 - I 1 )(R/^) r (0) n 2 s 

w i i a 11 1 x 

9 . 5V * 

+ (I x + I 2 " I 3 ) ^ 0 i + r (0) R()s " M 0. 

i i 1 


(16b) 


«RS v r.(L)(p.) a - -(ijR/i) qCL) p t - 777^ 777 + ^ 


(16c) 


^Rm r.(L)(s.) = (I 2 R /f) r (L) s + (I 3 - I 1 )(R/i) r . (L) fi 2 s 

W 1. 1 3. 11 -L-L 


+ di + i 2 - i 3 ) Q0 t + 


1 5v _ M * 

r . CLlR 5s . 0. 


(16d) 


where M^.., m| , M* and m|_ are the components of M" M|, M* and M=, 
1 1 i 1 

respectively. For the purpose of combining the treatment of the two 
variables v and w, we use the following matrix notations: 


< q i> 


s . 
1 


(q‘i) ; 



< u i> 



r i (e,x) (q ± ) 



<? ' 
H* 

O 

1 


<T 

H* 

O 

1 

(m i ) = 

0 “w. 

1 

(mj = 

1 

O 

1 


(p(e,x)) 


P v (e,x) 0 

0 P w ( £ ’ x 


( 17 ) 


12 


Substituting from equation (13) into equation (8) and making use 
of equation (17) results in 


00 



i=l 


+ m (X 


c) D r . (c,x) 
x i v ’ ' 


( q i> + r i (c,x) 


p(^»x) 


«i>a 


- 0. (18) 


By integrating the product of r n (c,x) and equation (18) from x = 0 to 
x = L and taking the limiting case that c approaches to zero, we obtain 


EIpL(q ) + m^ ) 
n n / , 

i=l 


f . (q . ) + pL(q ) 4- r^(0)(m ) (q ) + r^(L)(m ) (q ) = 0. 

m i H M n a n n / VM n / a n n' VM rra 


(19) 


Notice that use has been made of the orthogonality condition given by 
equation (11) and the integration formula given by equation (12). The 
constant f n ^ in equation (19) is defined as 


f . = 


P 1 

» -I - 




m 


cx) 


i. D^r. 


(x) + “ (x - c) D r . (x) y dx, 

m y v 1 I 3 

“ “ J 


( 20 ) 


which gives 


f u = f 15 = f,-, = f = f 

XX x^. _x nl n - :j 


= 0 f. ■ = 


PL 


f i = (-l) n ~ « B L (^ - 1) - 2a B L [ 1 - (-l) n ] l (n > 2) 
in m n n L n n u 


r- / \ ^ ^ Ld t / O _ 2 / n 

JL — ” -L ) Ul P -1^ V ^ ^ / J 

rn m L n n 


T - l T \ I O ^ 

v- i i-i j i i- 


A 

~J 


- 2Rd P|[l - (-l) n ] 1+ (-l^fVn L ' + 


13 


f . 

ni 


(-l) n+i — (L 


m 


2c) (a. p 

in 


a p.) - 4 r[i + (-1) 
n i 


n+i 


] (a p - a. p.) 

n n li ^4 


P 4 " P- 

' n 1 n 


( n , i § 2) 


PL j 5 

2ra [2 


(1 - 2c 


a / L) 
n 


(3 L 


(cqWl - 4) p 2 L 2 

n 3 7 K n 


+ R(a 8 l - 2) a p . 

v nn nn 


The last two terms of equation (19) are recognized as the dynamic coupl- 
ings given by equation (16). The partial derivatives in equation (16) 
are now determined from equation (5) by using equations (10) and (14). 
The results are written as follows: 


I djL 

R dp n 


mRP‘ J 


\ 

L 

i=l 


g ni P i 


I aV 
R 5p n 



8 ni P i 


I M 

R Bs 



e .s . 
6 ni l 


i=l 


I 

R 3s 


- - V 

= mRCT ^ 
i=l 


g -s . 
6 ni l 



0 


3v 



0 . 


(21) 


The formulations of g ni and g ni are shown in appendix B. 


We now have a set of four equations of motion, two of which are 
obtained from equation (19) by separating the variables and two of 
which are the third equation of equations (4a) and (4b). 


PL 


a (p / n 2 ) + 
n vr n 7 m 


[K(fyo 


oo 


1] 


+ I «„l'i ' V°> \ ' *„«■> 

i=l 


14 



b (s /n p ) 
n n ' 


t \ 4 ~ i 

— L) -r 

m K n 


I 3 - I x 


I 3 - 1 1 




0 r 2 (0) + 3 — r 2 (L) 

-.2 n - 2 n 


mi 


3 + ) C .S . 

n ni i 

i=l 


+ ar n (0) (0 n / H) + ar n (L)(0 n /fi) = r n (0) M“ + r n (L) M| (23a) 

n 


(P n /Q 2 ) + 50 n - br n (0)(s n /n) = M" 


(23b) 


( 0 n /fl 2 ) + 80 n - br n (L)(s n /fl) = , 


(23c) 


where 


a = pL/m + r 2 (0) I 3 /m £ 2 + r 2 (L) I 3 /m^ 
n n n ^ 


b n = pL/m + r 2 (0) I^/m ^ 2 + r 2 (L) Ip/mi0 


c . = f . + g . + (mR/mR)g . 
nx nx b nx v /5 nx 




a = (Ii + I 2 - I 3 ) / mRi & = (I 3 - I 2 )/Ix b = (1 + &) R 


.. / /? 
n; /j 


= (i 2 + I 2 - I 3 ) /mRi B = (I 3 - i 2 )/ii b = (1+5) R/i 


mRfPi M t 
n n 


✓w- 1 -y 

M/' — 1 — M 

0 mRn 2 ^ 0 

n n 


M" = — 2 M 
0 I-\Q 0 

n - n 


/S.u 1 * 

Mi = M- 

T p- f 

n mRn 2 iI 


m| = — - — m| m| 

0_ 2- X 0. 


M-' 

0 


n mRJTi n n I-Pr n 
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Notice that equations (22) and (23) represent the equations of motions 
of two inf inite-degrees-of -freedom systems of which the former is in the 
variable p n alone while the latter has three variables coupled together. 

In studying the dynamic stability or computing the frequencies of 
vibrations of a many-degrees-of-freedom system, it is often convenient 
to have the equations of motions in matrix forms. To this end, let us 
take a finite number of modes for each variable, say, N, we obtain 
readily from equations (22) and (23) 


•^2 Ap + Mp = E(t) 


(24) 


1 

s 

1 

s 


s 

F B 

— 1 

ra : na 
l 

+ - C 

a 

1 

«ra 

1 

+ K 

0 

0 


= F(t), 


(25) 


where 


P = 
N x 1 


Pi 


S 1 


0i 


01 

Ps 

s = 

s 2 

0 = 

02 

0 = 

02 

P N 

N x 1 

S N 

N x 1 

0 N 

N x 1 

0 N 


U = Identity Matrix D( ) = Diagonal Matrix 
N x N 


A = D(a ) Bj. = D (b ) R Q = D(r n (0)) 

NxN NxN NxN 


R L = D(r (L)) 
NxN 


M = 
NxN 


[c . ] 
ni 



[D(pV) - U] 



Bi 

0 

— 

0 


0 

aR o 

5r l 


Ki 

0 

0 

B 

0 

u 

0 

c 

-bR 

0 

0 

K 

0 

u 

0 

3N x 3N 

0 

0 

u 

3N x 3N 

o 

; br l 

0 

0 

3N x 3N 

0 

0 

u 
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Null Matrix 


K x = [c .] + 
ni 

N x N 


- In 


mr 


R 2 + 


- I- 




Rf + — KD(B 4 L 4 ) 
L m v n 


0 = 
N x N 


E (t) = M*R (e) + M-R (e) 
N x 1 * ° ^ 


F(t) 
3N x 1 


M*R (f) + M-R (f) 
0 O v ' 0 L v ' 

Mj(g) 

M|(g) 


In the expressions of E(t) and_F(t), we have introduced the load 
tion vectors, (e), (e), (f), (f ) , (g), and (g), such that 



r - 1 


— ■“ 


— 


H 

* 

s 


M* 

0i 


M* 

01 

M^e) = 

% 

M e (£) - 

M* 

> 

«>> - 

M* 

0p 


Hi 

N 


M* 


M* 

cA 

''N 

— 


M- 


M i 


M‘f 

01 

M| ( e) = 

M* 

# ^2 


M i 

S 

OQ 1 
II 

M* 

02 


M* 

% 


M- 

; u n 


M* 

0 N 


These vectors remain undetermined at present time. 


d is tr ibu- 


(26) 
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By transformation of variables, equation (25) can be reduced to the 
first order matrix differential equation [7], 


i y - Qy = G(t) , (27) 

where 



s 

0 





y = 

| 

Q = 

Q i Qz? 

u 0 

L j 

G(t) = 

0 

6N x 1 

S 

0 

0_ 

6N x 6N 

6N x 1 

F(t)_ 



0 -aB 2 R 0 -aB 2 R^ 


-B Ki 0 0 

Qi = "B~ 1 C = 

bR o 0 0 

Q 2 = -B _1 K = 

0 

1 

c 

o 

3N x 3N 

bR L 0 0 

3N x 3N 

o o -Eu 








B 2 = B" 1 = D(l/b n ). 


G. Solution of the Homogeneous Differential Equations 

Let the complementary solution of equations (24) and (27) be, 
respectively, 


p - PeV nt 

>4 

II 

(28a) 

y = Se Aflt . 


(28b) 
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Substitution of the assumed solution into the respective reduced differ- 
ential equations yields 

(M - 7 2 A)P = 0 (29a) 

(AU - Q)S = 0 (29b) 

and the characteristic equations 

|m - 7 2 a| = 0 (Nth degree polynomial of y 2 ) (30a) 

| AU - Q | = 0 (6Nth degree polynomial of A). (30b) 


Since the sub-matrices of matrix Q are either null or diagonal 
matrices except one, B^K.^, the 6Nth order determinant given by equation 
(30b) can be readily reduced to an Nth order determinant, 


| (A 2 + 5) (A 2 + SXA^i + K x ) + A 2 (A 2 + 8) abR 2 + A 2 (A 2 + S)abR 2 | = 0. 

(31) 

This gives a 3Nth degree polynomial of A 2 . For a symmetric system, all 
the barred symbols are equal to the unbarred symbols and R 2 = R^; the 
determinant of equation (13) degenerates to an Nth order determinant. 


(A 2 4- 5)(A 2 B 1 + K x ) + 2ab A 2 R 2 j = 0, 


(32) 


which is a 2Nth degree polynomial of A 2 . This can be seen from the 
relationship 


0 n = (-l) n 0 n n = 1,2, ... 


(33) 


obtained by eliminating the variable s n from equations (23b) and (23c). 
Hence, N variables of 0 n can be eliminated from the set of 3 n variables 
of a symmetric system. In equation (31) the constants a, a, b, b, 5, 5 
and the matrices B x , K x are directly related to the geometry and moments 
of inertia of the system and bending stiffness of the cable. The matrices 
R 2 , r£ and D(f3^L 4 ) can be written out immediately: 
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1 


1 

d 2 c 2 


d 2 (L - c) 2 

4 


4 

4 


4 

. 


• 

4 

_J 


*4 


D(f3jL 4 ) 


0 

0 


500.6 

3803 

14618 

39944 

89135 


Pn L 


1 

2 


(2n - 3) n 


n > 7 


R 2 




for a symmetric system . 
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H. Dynamic Stability and Free Vibrations of the System 
1 . Dynamic Stability Criteria 

From the solution of the homogeneous differential equations, the 
dynamic stability of a system which is free of external moments and internal 
damping force can be readily given as 
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I 


( 1 ) 

4 

> 0 

i = 1,2, ..., N for motions in 
rotation 

the plane of 

( 2 ) 

A 2 

< 0 

i = 1, 2, ..., 3N for motions 

out of the 


plane of rotation. 


2 . Free Vibrations of the System 

If a system is dynamically stable, its disturbed motions are 
vibratory and can be written in the forms: 


N 

i 

p = / (k cos y fit + k sin 7 fit) P 
' m n ns 'n n 

n=l 


(34) 


A fit- 


3n 

y ‘ Z Re { S n e ” 

n=l 


(A f s are positive, pure 


(35) 


imag inary) 

where Re { } means the real part of { }. In the above equations, P n and 
S n denote the eigen-vectors which are the solution vectors of 


(M - 7 2 A) P 
v 7 n n 


(A U - Q) S 
v n x n 


= 0 


n = 1 , 2, . . . , N 

(36) 

n = 1 , 2 , 3n, 

(37) 


respectively. The arbitrary constants k ni and k n2 ^ n cl the arbitrary 
magnitudes of the eigen-vectors S n can be determined for a given set of 
initial deformation and velocity of the cable as functions of x by using 
eigen-function expansion method. However, this is of little practical 
interest to us. 
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I. Special Case - Free Vibrations of Two Identical Stations 
Connected by a Massless, Flexible Cable 


Before we proceed to determine the responses of the system to external 
moments, let us apply the results of the previous section to a simple case. 
If the connecting cable is massless and flexible, it remains straight under 
centrifugal force. This is to say that the disturbed motions of the space 
stations produce only the two rigid-body modes of the cable. For free 
vibrations of a symmetric system connected by a massless and flexible 
cable, we substitute the following in equation (32), 



1 0 

1^ 

2 0 

Ai+ Is - I 3 \ 2 

2 0 

R 2 = Rf = 


B x = Tl 


, 2abR 2 = ( ) 


o L 

0 3 

x !r 

0 6 

o \ £ J 

0 6 


1 

2J X 0 

I 3 - I ! 

2 0 


II 


+ £ z 


, Ji = mRi Jp = 2mR 2 i/L 

/j 

0 6J 2 


0 6 
— _i 



and obtain the characteristic equation of the motions out of the plane of 
rotation 


where 


2fi(A) 

0 


0 

6f 2 (A) 


= 0 , 


(38) 


f n (A) = IjlsK 4 + [I 3 (I 3 " I 2 - Ii) + 2I 1 I 2 + I iJ n ] A 2 + (I 3 - I 2 )(I 3 " Ii + J n ), 

n = 1, 2. 

This gives two values of A for each mode from the equations 


f x (A) = 0 and f 2 (A) = 0. 
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Inserting 


I 

a * # 


2 

0 

1 

2Ji 

0 

0 

6 


0 

6J 2 


into equation (30a) yields 


7 n = ^V l3 ’ n = 1 , 2, 

for the motions in the plane of rotation. 


J. Responses of the System to External Moments 


Let us consider a system which is dynamically stable; i.e., the 
eigenvalues / n (n = 1, 2, . . . , N) are real and positive and the eigen- 
values Ah (n = 1, 2, . 3N) are pure imaginary. We now define the 

following: 

P = (PiP^ ... P^) modal matrix of matrix (A^M) or a row matrix of 
N x N the eigen-vectors of equation (29a) 

S = modal matrix of Q 
6N x 6N 


T 2 = D (7 2 ) = P -1 (A“ 1 M)P 
N x N n 


(39) 


A 

6N x 6N 


3 An 
0 



S _1 QS 


A x = D(A n ) 


sin FQt = D (s in y ut) 
N x N 11 


Afit 

e 

6N x 6N 


D (e 


Aq $7t > 




D(e '' n “ w ) 
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By making the transformation of variables 


p = Pq 


and 


y = Sz, 


( 40 ) 


(41) 


in equations (24) and (27), respectively, we obtain 


q + («r) 2 q = P“ 1 A“ 1 E (t) 


(42) 


and 


z - (flA)z = nS - 1 G(t). 


(43) 


It follows immediately that 


p = Pq = PT"’ 1 J sin m(t - t) P* 1 A“ 1 E(t) dx 
o 


(44) 


y = Sz = S / S -1 G(x) dr. 


(45) 


In the above, formulations of the solutions E(t) and G(t), as mentioned 
earlier, are yet to be determined. 

Before we proceed to determine the moment distribution vectors defined 
in equation (26), some considerations about the applied moments are neces- 
sary. To confine ourselves to the assumptions made earlier for the deriva- 
tion of the equations of motion, we require that (1) the resultant moment 
in the plane of rotation is small such that the change of the rotational 
velocity Q is negligible and (2) the resultant moment normal to the plane 
of rotation is also small so that the motions of the station normal to 
the orbital plane can be disregarded. Based on the above considerations 
and on the assumption that the moment of inertia of the system (mR 2 + mR 2 ) 
is much greater than the principal moments of inertia of the stations, we 
may assume that the disturbed motions of the stations due to the applied 
moments do not have appreciable effect on the motions of the whole system. 
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In other words, the applied moments which are used for the purposes of 
navigation, changing orbital plane, spinning up or down of the system, 
etc., are excluded. Although there is no means by which to determine the 
exact components of the applied moments on each generalized coordinate, 
it seems plausible to determine the loads distribution in a static equili- 
brium state and use it as an approximation in the solution of dynamic 
responses . 

Static Loads Distribution 


Let ^ and \|r be the angular displacements of the stations and u(x) be 
the deflection of the cable produced by the application of a unit moment 
M^. The work of the unit moment is equal to the sum of the change of 
potential energy of the stations and the strain energy of the cable. 


1 • \|r = AV + AV + £D. 


(46) 


Notice that the kinetic energy of the system is omitted under the assump- 
tion that the system is in static equilibrium configuration. We now 
assume that 


N 


u(x) 



c - r . (x) , 
1 1 


(47) 


of which the constants c^'s are to be determined from the condition 


€ = \|r - (AV + AV + AJ), 


(48) 


is a minimum; i.e., 


gr- = o n = 1, 2, ..., N. (49) 

n 

Some of the terms in equations (48) and (49) can be written using results 
already obtained: 
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cAV 2 \ 

7 — = ruR Q > g ,c. 
oc / , m 1 

n 1=1 


5AV -"2 ^2 

7 — = mR^n 
oc 

n 


g .c . 
m i 


( 50 ) 


N 


* = 


~ c r ( 0 ) 
£ n n v ' 


_ R , m 

aT - 1 r ~ (0) 


n=l 


(51) 


ALJ = 


1 

2EI 


^ . p c5 M 

„2 j 5 au 1 / A , x , 

M dx "t = — / M -r dx 

x dc El J x oc 

n n 

o 


77 A-i -c., 

El [_ , ni x’ 
i=l 


(52) 


where 


x 

M x =^ c n / pft 2 (c 
n=l o 


- i)[r n (x) - r n (|) ] d£ - mRn ? [r n (x) - r n ( 0 )], 


as given by equation (A-l). It is tedious to perform the integration given 
in equation (52); however, with the aid of reference 5 a closed form of h n -j_ 
can be obtained. This work is omitted here. 

Substituting from equations (50) to (52) into (49) results in a system 
of equations in c-^ 

(mR 2 Q 2 g . + mRRH 2 g . + -Jr h . )c. = 7 r (0) n = 1, 2, . N (53) 
x b m b ni El nr 7 1 £ n v 5 5 ' 


from which can be solved. We now replace the static deflection r n (x) 
by r n (x) P n (t) to represent the dynamic deflection, and rewrite the work 
expression as follows: 


N N 

M ? (t) ' M * Z I C n r „ (0) P „ <C) = Z M t/n< t >' 

n=l n=l 


K. 

Since ^ ( t) = — r n (0) P n (t), the above expression gives 




c . 

n 
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(54a) 


Similarly, 


Vv 




c 


n 


(54b) 


=M; C 

t, V n 


r n 

Jlr (L) 


M| = M; 5 
e n e n 


(54c) 


where c n = — 


n - Jr (0)" Cn * relationship is obtained by comparing equation 


(53) with the equation which generates c , 

<»R^g ni + fiR R fl=I n . + ± h nl ) c ± - | r n (L). 


ni 


(53a) 


However, the proposed static approach cannot be applied to the pitch moments, 
since the rotsional vibration of the cable is not considered. Therefore, the 
response to these moments cannot be treated at the present time. 


We now return to equations (25) and (26) and rewrite the load func- 
tions in the forms; 


E(t) = (M*(t) R o +M|(t) R l ) (c n ) (55) 

N x 1 v v 


F(t) 
3N x 1 


(M; ( t)R 0 +M|(t)R L )(c n ) 

0 

0 


where 



V 


(56) 


K. Discussion and Recommendations 

A great deal of computer time can be saved in computing equation (45) 
by using partitioned matrices to take advantage of the fact that the eigen- 
values and eigen-vectors of matrix Q are in conjugate pairs. For this 
purpose, let us partition the matrices as follows: 


S 

6N x 6N 


"" 1 

j QEA 1 , 

- j QEAj_ 

w 7 

_ 2 , 

E 


L 

3N x 3N 


a 

M- 

v 


-l 

S 

6N x 6N 


IIjJ JL 

•V I ^ 

111 , n 


n = (§, t)» O 

3N x 3N 


(57) 


(58) 
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where M ~" above a letter denotes its conjugate; a, y are N x 3 n sub- 
matrices; and T), £ are 3N x N sub-matrices. It is easy to prove that 
matrix II can be computed directly from the equation 


n = (ZAj_ + ZAi Z 1 Z )“ 1 LA X Z " 1 


3N x 3N 


(59) 


Now, equation (45) becomes 


s = 2Re {crH(t) } s = 2Re {jaA;,. H(t) } 

N x 1 N x 1 


N x 1 


= 2Re (pH(t) 


0 = 2Re {j^ H(t) } 

N x 1 


(60) 


0 = 2Re (vH(t) } 

N x 1 


0 = 2Re {jvAj. H(t) ), 

N x 1 


where Re denotes the "real part of" and 


H(t) = I e jAia(t ' T) |(M*(t)R o + m|(t)R 0 ) (c n ) d- 


(61) 


From equations (44) and (52), 


p = PT" 1 / sin ar(t - t)P - 1 A -1 (M', (t)R + M*(t)R ) (c ) dx (62) 

. / V o iLr u n 


Hr 


and 


p = fiP / COS flr(t - t) P _ 1 A _1 (M“(t)R 0 + M^.(t)R o ) (c n ) dr. 


28 



Finally, the Euler angles and their time derivatives of the space sta- 
tions are 


r-d ^ 
II 

<v 

^ 

11 

o 

i R s 
o 

® = Z0 


II 

^ 1 


CD 1 

II 

bo ill-* 

iR L S 

0 = Z0 


-€*-• 

II 

^ I*- 1 


CD* 

II 

bo 

<v 

• • 

0 = Z0 

(63) 

* 1 

* = - 
£ 

ZR L P 

0 " - 

I 

<r l s 

$ = i0 



where i=(ll...l),alxn row matrix. 

In conclusion, the following problems are suggested for further 
s tudy : 


(1) Refine the approach used for the determination of the 
load distribution factor . 

(2) Obtain solutions on other types of loadings such as reel- 
out and reel-in of the cable, spin-up and spin-down of 
the stations, movements of the astronauts, etc. 

(3) Take into consideration the torsional stiffness of the 
cable and other end conditions of the cable. 

(4) Study other configurations of the stations, such as two 
stations connected to a massive hub. 
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APPENDIX A 


Determination of M(x,t) and N(x,t) and Their Derivatives 


As observed from Figure A-l, the bending moment at x due to the 
centrifugal force is 


M(x,t) = - J pn 2 (c - g)[u(x,t) - u(|,t)] d | - mRH 2 [u(x , t ) - u(0,t)] x< c 
o 


M(x,t) = - J pfi 2 (g - c)[u(x,t) - u(|,t)] d| - mRP 2 [u (x, t) 
x 

Since x = c is the center of rotation, 


u(L, t) ] c< x <L. 
(A-l) 


mR + pc 2 = mR + « p(L - c) 2 . 


(A- 2) 


It can be proved readily that 





2cx) 


S 2 u 

Bx 2 


+ 2(x 



- mRp 2 


B 2 u 

Sx 2 5 


0 < x < L. 


(A- 3) 


The axial force in the cable is 


,<> - r 


N(x,t) = / p£T (c - |) dg + mRfr, x < c 


N(x,t) = J pfi 2 (i - c) dg + mRQ 2 , x > c. 
x 


(A- 4) 
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APPENDIX B 


The Potential Energy 


Substituting equations (13), (14) and (15) into equation (5) results 


rn 


N N 


-£R(q) = 1 
R 2 2 



f Gr) Gr) + 7 r n (0) r l (0) 


n, i=l o 


q q . 


N N 

- \ d^cqj - d Y [r .<C) - 2] , 2 , 1 + | £ [.„(.) + 12a n P n ] q = 
i=3 n=3 


N N 


+ - 



9 , / [a .(c) + 4(p 4 a.p. - pfcc 6 )](P 4 - Pf)" 1 q q. 

2 /-j Z_j ni v i K i ir ir VK n 

n, i=3 


N N 


4 



\ r (0) r . (0) q q . . 
i n v i v H n^i 


n, i=l 


(B-l) 


, I 

R 2 2 


N N L 

P /dr \ /dr ,\ . 

/Gr; Gr) d * + 

n, i=l c * 



q q * 


N 

= \ d 2 (L - c)q 2 + d V [(-1) L 2 + r (c)] q^ 


2 I_A 

i=l 


N 


+ i ) [4a (3 (3a (3 + 1) - a (c)] q 2 

8 / , n n n n 


nn n 
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(B-2) 


+ “ ) > [(-1) 4(p 4 a.p. - pfep) - a .(c)] q q. 

2 / , / , VK n i K i K i n^n' m 

n, i=3 


N N 


+ I 



- r (L) r . (L) q q . , 

- n i v 7 n n M i 


n, i=l 


where 


dr. dr d 2 r. d :j r d 2 r d 3 r. 


dT ' Pt r i 3T 1 - d^“ + dP* " * 1 


dr Air 

a (x) = 3r — — + x (—*■ 
nn n dx \dx 



d'-r 


" 2xr n dP* ‘ C dF^ d^ + (X/P n> (aTt) 


Let us denote 


N 


N 


~7 = mRa 2 ) g .q. 

R c>q L ni 1 

i=l 


- ^r~ = mRa 2 ; g„ . q . , 


A n 


R oq 


ni l 


i=l 


where the formulas of g n ^ and g n ^ are tabulated as follows; 


^in ^ni 


gn= R/Jt gip = dcR/ £ 


g = 2R/£ (n > 2) 
in 


g 2 . 2 . ~ Rd 2 c(l + c/ f ) 


g 2n = Rd [2c/ i + 2 - r n ( c )] 


(n > 2) 


a . (c)+4(p 4 a.p. - pfa p ) , - 

m n i H i n 2 n 4 


mi 


„ /I „ A 




R g = y a (c)+3a p + 
6 nn | 4 nn nm 


^ | 4 > 



(n > 2) 


^in ®ni 


gn — R/ H §12 — “d (L-c)R/ £ g in “(“!) 2R/ i 


g 22 = Rd^(L-c) (1 +-:- N ) £ = Rd[(-l) n 2(L-c)/i + (-l) n 2 

i J 2n 


g ni 


40%. p. - pfcc p )(-l) 

K n i i i iv n 


n+i 


6 4 - p? 

n i 


a . (c) . 

— — + (-i) 4/i 


R 


g = [a p (3a p + 1) - 7 a (c) + 4/^] R. 
6 nn L n n n H n y 4 nn 


+ r n (c>] 
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